Complementary Constraints on Brane Cosmology 



J. S. Alcaniz* 

Departamento de Astronomia, Observatorio Nacional, 20921-400 Rio de Janeiro - RJ, Brasil 

Zong-Hong Zhu^ 

Department of Astronomy, Beijing Normal University, Beijing 100875, China 
(Dated: February 2, 2008) 

The acceleration of the expansion of the universe represents one of the major challenges to our 
current understanding of fundamental physics. In principle, to explain this phenomenon, at least two 
different routes may be followed: either adjusting the energy content of the Universe - by introducing 
a negative-pressure dark energy - or modifying gravity at very large scales - by introducing new 
spatial dimensions, an idea also required by unification theories. In the cosmological context, the 
role of such extra dimensions as the source of the dark pressure responsable for the acceleration 
of our Universe is translated into the so-called brane world (BW) cosmologies. Here we study 
complementary constraints on a particular class of BW scenarios in which the modification of gravity 
arises due to a gravitational leakage into extra dimensions. To this end, we use the most recent 
Chandra measurements of the X-ray gas mass fraction in galaxy clusters, the WMAP determinations 
of the baryon density parameter, measurements of the Hubble parameter from the HST, and the 
current supernova data. In agreement with other recent results, it is shown that these models provide 
a good description for these complementary data, although a closed scenario is always favored in 
the joint analysis. We emphasize that observational tests of BW scenarios constitute a natural 
verification of the role of possible extra dimensions in both fundamental physics and cosmology. 

PACS numbers: 98.80.Es; 04.50.-|-h 



I. INTRODUCTION 

The idea of a dark pressttre-dominated universe is usu- 
ally associated with unknown physical processes involv- 
ing either new fields in high energy physics or modifi- 
cations of gravity at very large scales. If one chooses 
to follow the former route, then the two favorite candi- 
dates for this mysterious component are the energy den- 
sity stored on the true vacuum state of all existing fields 
in the Universe, i.e., the cosmological constant (A), and 
the potential energy density associated with a dynamical 
scalar field (</)), usually called dark energy or quintessence 
(see, e.g., |1| for some recent reviews on this topic). 

The second route in turn is predominantly related to 
the existence of extra spatial dimensions, an idea that is 
required in various theories beyond the standard model of 
particle physics, especially in theories for unifying grav- 
ity and the other fundamental forces, such as superstring 
or M theories. In the eleven-dimensional supergravity 
model of Hoi'ava and Witten 0], for instance, the or- 
dinary matter fields are assumed to be confined in a 
submanifold (brane) immersed in a higher dimensional 
space, usually named bulk. An important development 
of this idea was subsequently given by Randall and Sun- 
drum Q who showed that, if our three-dimensional world 
is embedded in a four-dimensional anti-de Sitter bulk, 
gravitational excitations arc confined close to our sub- 
manifold, giving rise to the familiar law of gravity 
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(see also |j| for further discussions). 

From the observational viewpoint, the most important 
aspect asssociated with these scenarios resides on the 
fact that some of their versions can lead to a late time 
acceleratinguniverse, in agreement with supernovae ob- 
servations For this reason (and also motivated by 
a possible explanation for the hierarchy problem from 
extra-dimension physics), braneworld (BW) cosmologies 
has become a topic of much interest recently. Sahni and 
Shatnov 0], for instance, proposed a new class of BW 
models which admit a wider range of possibilities for the 
dark pressure than do the usual dark energy scenarios. 
As shown in Ref. j^, for a subclass of the parameter 
values the acceleration of the Universe can be a transient 
phenomena, which could help reconcile the supernova ev- 
idence for an accelerating universe with the requirements 
of string/M-theory More recently, Maia et al. Q 
showed that the dynamics of a dark energy component 
parameterized by an equation of state p — cup (the so- 
called XCDM parametrization) can be fully described by 
the effect of the extrinsic curvature of a FRW universe 
embedded into a five-dimensional, constant curvature de- 
Sitter bulk. 

Another interesting scenario is the one proposed by 
Dvali et al. 0, which we will refer to it as DGP 
model. It describes a self-accelerating 5-dimensional BW 
model with a noncompact, infinite- volume extra dimen- 
sion whose dynamics of gravity is governed by a competi- 
tion between a 4-dimensional Ricci scalar term, induced 
on the brane, and an ordinary 5-dimensional Einstein- 
Hilbert action (see ilfl] for details). DGP models have 
been sucessfully tested in many of their observational 
predictions, ranging from local gravity to cosmological 
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observations [11111111111111 (see, however, [HITl). 
From the theoretical viewpoint, the consistency of their 
self-accelerating solutions is still a matter of debate in 
the current literature (see, e.g., [13 )■ 

This paper aims at placing new observational con- 
straints on DGP models from the current X-ray observa- 
tions of rich clusters of galaxies and type la supernovae 
(SNe la) data. We emphasize that this particular com- 
bination of observational data constitutes an interesting 
and complementary probe for testing the viability of cos- 
mological scenarios because while X-ray data are very 
effective to place limits on the clustered matter (dark 
matter) the new SNe la sample tightly constrains the un- 
clustered component (dark pressure). To this end, we use 
the latest Chandra measurements of the X-ray gas mass 
fraction in 26 galaxy clusters, as provided by Allen et 
al. ^3 along with the most recent determinations of the 
baryon density parameter, as given by the WMAP team 
[l^. the latest measurements of the Hubble parameter, 
provided by the HST key project [53, and the so-called 
old set of 157 SNe la, recently published by Riess et al. 
The analysis performed here, therefore, updates the 
results of Ref. 1211 a nd complements the SNe la study 
presented in Ref. |2^. 



II. DGP MODELS: BASIC EXPRESSIONS 

In DGP models, the presence of an infinite- volume ex- 
tra dimension modifies the Friedmann equation in the 
following way JXJj 
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where p is the energy density of the cosmic fluid, k = 

0, ±1 is the spatial curvature, Mpi is the Planck mass 
and Tc = Mpj/2M| (M5 is the 5-dimensional reduced 
Planck mass) is the crossover scale defining the gravita- 
tional interaction among particles located on the brane, 

1. e., for scales r < rc, the gravitational force experienced 
by two punctual sources is the usual 4-dimensional l/r^ 
force whereas for distance scales r > rc the gravitational 
force follows the 5-dimensional 1/r^ behavior. From the 
above equation we find that the normalization condition 

is now given by flk + [a/ fire + y/^r^ + ^^m] = 1, where 
ri,n and Qk are, respectively, the matter and curvature 
density parameters (defined in the usual way) and 
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is the density parameter associated with the crossover 
radius rc- For a flat universe, the normalization condi- 
tion above reduces to fii-^ = (1 — flm)^ /4. As noticed 
in Ref. 0, the above described cosmology can be ex- 
actly reproduced by the standard one plus an additional 
dark energy component with a time-dependent equation 
of state parameter Lu^^f{z) = 1/Q{z, il^, fti-J — 1, where 



-3 -I- 1) and x' = {1 + z)"i (see l23| for a 
discussion on time-dependent parametrizations for u). 

To perform our statistical analysis in the next section 
two observational quantities are of fundamental impor- 
tance, namely, the angular diameter [dA(-z)] and luminos- 
ity distances [cZl(-z)] - intrinsically related, in a homoge- 
neous and isotropic universe, by dA{z){l + z)^ = di^{z). 
From the above equations, it is straightforward to show 
that 



jDGP 
"A 



iz) 



H- 



{i + z)\nk\^/^ 



(3) 



xSk 



dx 



X^J^iflj 



where the function Sk is defined by one of the follow- 
ing forms: Sk{r) = sinh(r), r, and sin(r), respectively, 
for open, flat and closed geometries. The dimcnsionless 
function x) is given by 

(4) 

where j stands for m, rc and k. As one may check, for 
rife = 0, the limit 1/rc {fl^^ — *■ 0) provides 
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which is the standard (SCDM) prediction for the angular 
diameter distance dAiz). 



III. CONSTRAINTS ON DGP MODELS 



The X-ray gas mass fraction test [/gas(-z)] was first in- 
troduced in Ref . l23 a nd further developed in Ref. |25| 
(see also [2ll[2l[23l23). This is based on the assumption 
that rich clusters of galaxies are large enough to provide a 
fair representation of the baryon and dark matter distri- 
butions in the Universe |2^. Following this assumption, 
the matter content of the Universe can be expressed as 
the ratio between the baryonic content and the gas mass 
fraction, i.e., f2ni oc rib//gas- Moreover, as shown by 



Sasaki |24| since /gas 
defined as [lal 
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the model function can be 
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where the bias factor 6 is a parameter motivated by gas 
dynamical simulations that takes into account the fact 
that the baryon fraction in clusters is slightly depressed 
with respect to the Universe as a whole 30], the term 
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FIG. 1: The model function /™° as a function of the redshift 
for some selected values of and Q.th'^ = 0.0224, b = 0.824 
and h — 0.72. The value of Qm is fixed at 0.3, as suggested 
by dynamical estimates |3lll . The current concordance model, 
i.e., a fiat scenario with 70% of the critical energy density 
dominated by a cosmological constant, is also shown for com- 
parison. 



0.19Vh stands for the optically luminous galaxy mass in 
the cluster and the ratio dj^^^ (zi) / cff^'^^ (zi) accounts 
for deviations in the geometry of the universe (here mod- 
elled by the DGP model) from the default cosmology used 
in the observations, i.e., the SCDM model (see O 113 
for more observational details). 

In order to constrain the cosmological parameters flm 
and r^r^ we use the latest Chandra measurements of the 
X-ray gas mass fraction in 26 dynamically relaxed galaxy 
clusters (0.07 < z < 0.9), as provided by Allen et al. [l8| . 
The confidence intervals are determined by using the con- 
ventional minimization with the Gaussian priors on 
the baryon density parameter, Qhh'^ = 0.0224 ± 0.0009 
[ill, on the Hubble parameter, h = 0.72 ± 0.08 |23|, and 
on the bias factor, b = 0.824 ± 0.089 Hi, i.e.. 
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In the above expression, f^'^{z\) is given by Eq. (6) and 
/gas, i is the observed values of the X-ray gas mass fraction 
with errors j. To plot two-dimensional confidence 
contours, we have projected our five-dimensional parame- 
ter space (defined by the vector p = {fi™, h, ^bh^, b}) 
into the plane flm — ^r^, which is similar to marginal- 
ize over the parameters h, fl^h^ and b by defin- 



ing the probability distribution function £(i7,„, i7i.^) = 
J dhdinbh^)dbe-^' (see, e.g., H). 

Figure 1 shows the behavior of f^°'^ as a function of 
the redshift for some selected values of fti-^ and flbh^ — 
0.0224, b = 0.824 and h ^ 0.72. The value of 17™ is fixed 
at 0.3, as suggested by dynamical estimates on scales up 
to about 2h~^ Mpc 31]. For the sake of comparison, the 
current favored cosmological model, namely, a flat sce- 
nario with 70% of the critical energy density dominated 
by a cosmological constant (ACDM), is also shown. We 
note from this figure that with a larger sample of high-z 
X-ray measurements {z > 1) it will be possible to com- 
pletely distinguish between DGP and ACDM scenarios. 
In Fig. 2a we present the first results of our statisti- 
cal analysis. There, it is shown the confidence regions 
(68.3%, 95.4% and 99.7%) in the plane - by con- 
sidering the 26 X-ray data points described above. The 
best-fit parameters for this analysis are flm = 0.29 and 
firo — 0.31 with x^/t^ — 1-07 is defined as degrees 
of freedom). In particular, this relative value of is 
slightly larger than the one found in Ref. 0| for the 
ACDM case with arbitrary curvature, i.e., x^/i^ — 1.02. 
At 95.4% confidence limit (c.l.) we obtain the intervals 
0.23 < 0„ < 0.37 and 0.25 < Q,, < 0.37 whereas by 
restricting our analysis to the flat case, we note that the 
data favour a lower value of the matter density param- 
eter, i.e., ~ 0.23 {Qr, ~ 0.148) with xLJ'' - 1-2- 
Such a value, however, is inside the la interval of the 
WMAP and other recent estimates of the quantity of 
matter in the Universe, fim = 0.27 ±0.04 j^, and there- 
fore does not imply any possible conflit between the pre- 
dictions of the model and the independent measurements 
of (see for a discussion on this point). 



B. Complementary Constraints 

By comparing our Fig. 2a with Fig. 2 of Ref. p3 | 
(also shown here as Fig. 2b for the sake of comparison), 
it is clear that X-ray and SNe la data provide orthogonal 
statistics in the plane flm — f^r,,- This, therefore, sug- 
gests that possible degenenacies between these parame- 
ters may be broken by combining these two data sets in a 
joint statistical analysis. To perform such an analysis we 
make use of the so-called gold SNe la sample of Riess et 
al. H, which consists of 157 events distributed over the 
redshift interval 0.01 ^ z 1.7. This particular sample 
constitutes the selection of the best observations made 
so far by the two supernova search teams along with 16 
new events observed by the HST. Our app roach for SNe 
la data is based on Ref. (see also |32|'). 

The results of our joint analysis are shown in Figure 
2c. Note that the combination of these data sets leads to 
tight constraints on the — f^ic plane. A comparison 
with Fig. 3 of Ref. :21] (which used the then available 
172 SNe la taken from Tonry et al. [13 plus 9 X-ray clus- 
ters from Allen et al. shows that the current anal- 
ysis reduces considerably the area corresponding to the 
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FIG. 2: The results of our statistical analyses: Confidence regions (68.3%, 95.4% and 99.7%) in the fim — ^^rc plane by 
considering: a) the latest Chandra measurements of the X-ray gas mass fraction in 26 galaxy clusters (0.07 < z < 0.9) plus 
determinations of the baryon density parameter and measurements of the Hubble parameter; b) the so-called Gold sample 
of Riess et al. - taken from [2^ and shown here for the sake of comparison. Note that X-ray and SNe la data provide 
orthogonal statistics in the plane Q.m — f^rc; c) Joint X-ray -|- SNe la + Q.iyh'^ + Ho analysis. This combination of data provides 
Q,-m = 0.34jlQ Qg and f2r^ = 0.26 ± 0.04 (at 95.4% c.l.) and clearly favors a closed universe. 



TABLE I: Recent estimates of the crossover radius 
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confidence intervals. As expected, the parameter space 
now is reduced relative to our former analysis, with the 
best-fit scenario occurring at Vl^ = 0.34 and Vt-^^ = 0.26 
{x^ /v = 1.11). These values correspond to an accelerat- 
ing universe with deceleration parameter Qo — —0.8 and 
a total expanding age of to — 10.1h~^ Gyr. At 95.4% 
c.l. we found the following intervals: 0.29 < ftm S 0-41 
and 0.22 < il^^ < 0.30. Note that this 2a interval for 



leads to an estimate of the crossover scale rc in terms of 
the present Hubble radius H^^ [see Eq. (2)], i.e., 

= 0.98t^:;5?i/-\ (8) 

which is slightly smaller than (but in agreement with) 
the value obtained in Ref. = imto fe^-^ (99% 

C.L). If we restrict our analysis to a flat geometry, i.e., by 
imposing the normalization condition f2r^ — (1 — r2m)^/4, 
we obtain flm — 0.22 {fi^^ ~ 0.152), which corresponds 
to a 9.6/i~^-Gyr-old universe with Qo = —0.45 and Tc — 
1.28H-\ 

At this point we compare our estimates of the crossover 
scale Tc with other recent determinations of this quantity 
from independent methods. We note a reasonable agree- 
ment among them. For instance, in Ref. (13 Deffayet et 
al. used the then available SNe la -I- CMB data to find 
Tc ~ 1.4_ff~^ for a flat model with = 0.3, in agree- 
ment with Avelino and Martins who used a larger 
sample of 92 SNe la events. The current measurements 
of the angular size of high-z sources require an accel- 
erating universe with a crossover radius of the order of 
Tc ~ 0.947?~^ [iJI whereas the statistics of gravitation- 
ally lensed quasars implies Vc ~ 1.76H~^ UM- The less 
concordant (but not in desagreement) among the current 
estimates for Vc comes from age measurements of high-z 
galaxies, which requires Vc < 2.04:H~^ ^M- These re- 
sults, along with the main estimates of the present paper, 
are summarized in Table I. 
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IV. FINAL REMARKS 

There is increasing evidence for an accelerating uni- 
verse from various astronomical observations. However, 
understanding the acceleration mechanism based on fun- 
damental particle physics is still one of the most impor- 
tant challenges in modern cosmology. An unknown dark 
energy component with negative pressure has usually 
been invoked as the most feasible mechanism for the ac- 
celeration although effects arising from extra dimension 
physics can mimic dark energy through a modification on 
the Friedmann equation |4j .6, .8 , 34] . In this paper, we 
have focused our attention on a specific self-accelerating 
five-dimcnsional braneworld scenario, the so-called DGP 
model 9J. 

We analyzed the scenario by using the most recent 
Chandra measurements of the X-ray gas mass fraction in 
galaxy clusters, the WMAP determinations of the baryon 
density parameter, measurements of the Hubble param- 
eter from the HST, and the current supernova data. As 
shown, the model provides a good description for the 
so-called gold SNela sample of Ref. |^. However, in 
order to explain the X-ray data of clusters in the frame- 
work of DGP model, a closed universe is necessary (a 
similar result is also obtained by fitting SNela data to 
the standard ACDM model. In this case, CMB data is 
generally used to match a flat universe.). The combi- 



nation of the above mentioned data sets leads to very 
tight constraints on the flm — plane. At 95.4% c.l. 
we found the following intervals: 0.29 < flm ^ 0.41 and 
0.22 < r^r^ ^ 0.30, which gives a closed universe with the 
curvature of —0.39 ^ ilfc ^ —0.18. It is worth mention- 
ing that we might make heavy use of the X-ray gas mass 
fraction in clusters, which further prefers to a closed uni- 
verse in DGP model. This kind of analysis depends on 
the assumption that the /gas values should be invariant 
with redshift, which has been criticised by a minority 
of workers in the field. For example, a recent compar- 
ison of distant clusters observed by XMM-Newton and 
Chandra satellites with available local cluster samples 
seems to indicate a possible evolution of the M-T re- 
lation with redshift, which may be indicating that the 
standard paradigm on cluster gas physics needs to be re- 
vised [s^. Therefore, to pin down the DGP model from 
this kind of observations, a more detailed study on X-ray 
gas mass fraction in galaxy clusters is necessary. 
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